home *** CD-ROM | disk | FTP | other *** search
/ Aminet 1 (Walnut Creek) / Aminet - June 1993 [Walnut Creek].iso / aminet / misc / sci / ephem_src_4_28.lha / obliq.c < prev    next >
C/C++ Source or Header  |  1992-04-17  |  421b  |  22 lines
  1. #include <stdio.h>
  2. #include "astro.h"
  3.  
  4. /* given the modified Julian date, mjd, find the obliquity of the
  5.  * ecliptic, *eps, in radians.
  6.  */
  7. obliquity (mjd, eps)
  8. double mjd;
  9. double *eps;
  10. {
  11.     static double lastmjd = -10000, lasteps;
  12.  
  13.     if (mjd != lastmjd) {
  14.         double t;
  15.         t = mjd/36525.;
  16.         lasteps = degrad(2.345229444E1
  17.             - ((((-1.81E-3*t)+5.9E-3)*t+4.6845E1)*t)/3600.0);
  18.         lastmjd = mjd;
  19.     }
  20.     *eps = lasteps;
  21. }
  22.